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We investigate the behavior of oxygen vacancies in three different metal-oxide semiconductors (rutile and 
anatase Ti 02 , monoclinic WO 3 , and tetragonal Zr 02 ) using a recently proposed hybrid density-functional 
method in which the fraction of exact exchange is material-dependent but obtained ah initio in a self-consistent 
scheme. In particular, we calculate charge-transition levels relative to the oxygen-vacancy defect and compare 
computed optical and thermal excitation/emission energies with the available experimental results, shedding 
light on the underlying excitation mechanisms and related materials properties. We find that this novel 
approach is able to reproduce not only ground-state properties and band structures of perfect bulk oxide 
materials, but also provides results consistent with the optical and electrical behavior observed in the corre¬ 
sponding substoichiometric defective systems. 

I. INTRODUCTION 

Oxygen vacancies represent a prototypical intrinsic 
point defect in metal-oxide materials, where they are in¬ 
evitably formed throughout the synthesis process, as a 
result of the exposition of the samples to high tempera¬ 
tures and strongly reducing atmospheres. 

The presence of oxygen vacancies substantially al¬ 
ters the optical and electrical properties of metal-oxide 
semiconductorsand a correct interpretation of exper¬ 
iments requires these modifications to be properly taken 
into account. Moreover, defects control and engineering 
play a fundamental role in a wide range of applications, 
such as in catalysis and photocatalysis,— information 
technology,— gas sensors^ and smart windows.— 

The standard approach to point defects computational 
modeling in materials is rooted within density-functional 
theory (DFT) and the supercell approximation.— The 
simplest interpretation of the spectroscopic features of 
defective materials may be provided by analysis of the 
corresponding Kohn-Sham (KS) electronic band struc¬ 
tures. However, it is well known that KS energy eigen¬ 
values represent a very crude approximation to quasi¬ 
particle energies probed in spectroscopies. The limits 
of such approximation noticeably manifest themselves in 
the failure of standard DFT in predicting accurate elec¬ 
tronic band gaps. An exact xc functional, providing ex¬ 
act total energies, would be able to yield exact defect 
states through comparison of defect formation energies 
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corresponding to variations of the charge state of the de¬ 
fect {charge-transition levels, CTLs)— In principle, 
CTLs can be rigorously put in correspondence with ex¬ 
perimentally measurable excitation and emission energies 
in defective semiconductors and insulators in an anal¬ 
ogous spirit, the delta-self-consistent-field method allows 
one to compute electronic excitations in finite systems 
from total-energy differences.— In practice, of course, the 
error due to an approximate form of the xc term would 
affect total energies and hence the computed CTLs. No¬ 
tice in fact that standard local/semilocal DFT function¬ 
als often predict defect states to be resonant with the 
bulk band edges, even when experiments indicate a deep 
donor or acceptor behavior; in this case, computed total 
energies, and the corresponding CTLs, would be largely 
incorrect— Hence, improving the functional with respect 
to local/semilocal ones is a prerequisite for the applica¬ 
tion of the method. Promising results have been achieved 
by adopting exact-exchange hybrid DFT functionals^^ 
for the description of exchange and correlation— Such 
functionals have been proven to yield band gaps in more 
quantitative agreement with experiments, as well as to 
help correcting for the self-interaction error inherent in 
local/semilocal approximations. 

In the frame of hybrid DFT calculations on extended 
systems, a connection between the fraction of exact ex¬ 
change to be admixed in the hybrid functional and the 
dielectric constant of the material has been proposed— 
The obtained functional has proved capable of predict¬ 
ing band gaps, structural properties and total ener¬ 
gies for semiconductors in quantitative agreement with 
experiments— In particular, in our recent worl— we 
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tested the performance of this approach in calculating 
various properties of metal-oxide semiconductors. 

Motivated by this success, we devote the present work 
to an assessment of the efficiency of our method in 
describing defective oxide semiconductors, with partic¬ 
ular concern to oxygen vacancies, which constitute a 
prototypical case study in this kind of materials. In 
this respect, we investigate the properties of three dif¬ 
ferent oxide semiconductors of great both fundamental 
and applicative interest: rutile and anatase titanium 
dioxide (Ti02), tungsten trioxide (WO 3 ) in its room- 
temperature (RT) 7-monoclinic crystal structure, and 
tetragonal zirconium dioxide (Zr02). 

The paper is organized as follows. In Section M de¬ 
tails on the performed DFT calculations are given, the 
dielectric-dependent hybrid DFT method is discussed 
and the charge-transition levels formalism is reviewed. In 
Section uni results are discussed and interpreted in the 
light of the available experiments on reduced oxide sam¬ 
ples. Finally, the most important conclusions are sum¬ 
marized in Section HVl 


II. COMPUTATIONAL APPROACH 

A. DFT calculations 

All the calculations were performed using the CRYS- 
tal09 code^Sii^i based on linear combination of atomic 
orbitals and Gaussian-type basis sets. Adoption of small¬ 
sized localized atomic-like functions for representing KS 
orbitals (compared, e.g., to much larger plane-wave ba¬ 
sis sets) allows one to considerably reduce the computa¬ 
tional cost of hybrid-functional calculations almost down 
to that of standard local/semilocal DFT calculations. 
This in turn allows to employ moderately large super¬ 
cells (containing roughly a hundred atoms) for modeling 
defective materials, requiring a reasonable computational 
effort also when performing structural optimizations. 

In the case of Ti 02 , calculations were performed in the 
all-electrons scheme, using the basis set from Ref. for 
Ti. Small-core effective-core pseudopotentials were em¬ 
ployed for modeling core electrons for W and Zr atoms2^ 
in WO 3 and Zr 02 , while valence electrons were described 
using the basis sets defined in RehM (5p, 5d, 6sp W elec¬ 
trons in the valence) and Ref. (dd, 5sp Zr electrons 
in the valence), respectively. The O atom was always 
treated at the all-electron level, using the basis set from 
Ref. m for Ti 02 and WO 3 , and from Ref. [13 for Zr 02 . 

We considered bulk 72-atom 2x2x3 and 96-atom 
2-\/2 X X 2 supercells for rutile and anatase Ti 02 , 
respectively. For RT 7-monoclinic WO 3 , we adopted a 
model supercell comprising 64 atoms, obtained by dou¬ 
bling the primitive cell along the a crystallographic axis. 
For tetragonal Zr 02 , a 108-atom 2x2x3 supercell was 
employed. The O vacancy was modeled by removing one 
O atom from the corresponding supercells. The Brillouin 
zone was sampled using 2x2x2 F-centered Monkhorst- 


Pack grids^ corresponding to 6 or 8 fc points in the 
irreducible wedge. 

Structural optimizations were performed by allowing 
all the atoms in the cell to relax their positions, keeping 
the lattice parameters fixed to that of the optimized bulk 
cell. Convergence thresholds in geometry optimizations 
were set at their standard values in the CRYSTAl09 code: 
the defined thresholds for the maximum and the root- 
mean-square of the energy gradients (atomic displace¬ 
ments) are 0.00045 a.u. (0.00180 a.u.) and 0.00030 a.u. 
( 0.00120 a.u.), respectively!^ 

B. Dielectric-dependent self-consistent hybrid functional 

The xc potential was approximated within the general¬ 
ized KS scheme^ using a full-range hybrid functionally 
in which the fraction of exact exchange a is computed 
self-consistently with the static electronic dielectric con¬ 
stant CooIn the framework of many-body pertur¬ 
bation theory (and in particular within the GW scheme), 
it can be proven that a ~ l/coo, provided that micro¬ 
scopic components of the dielectric function are averaged 
out over the cell and dynamical effects are neglected (in 
which case the GW self-energy reduces to the Hedin’s 
Coulomb-hole-plus-screened-exchange approximatiory^) . 

The electronic dielectric constant Coo was com¬ 
puted within first-order perturbation theory using the 
coupled-perturbed Kohn-Sham approach implemented in 
CRYSTAL09j^ allowing for evaluation of the xc contri¬ 
bution to electronic screening beyond the random-phase 
approximation. 

The dielectric-dependent self-consistent hybrid func¬ 
tional was defined for the bulk pristine crystal, and then 
applied to the corresponding defective system. This pro¬ 
cedure is justified by the facts that (i) the method proved 
able to provide electronic gaps, structural properties and 
ground-state energies in quantitative agreement with ex¬ 
periment for bulk systems^yy and (ii) the presence of 
point defects is expected to negligibly affect the dielec¬ 
tric properties in the dilute limit of defect concentration 
which is ultimately addressed. 

The self-consistent procedure yielded the following val¬ 
ues for the exchange fraction: 15.2% for rutile and 18.4% 
for anatase Ti 02 , 21.9% for RT monoclinic WO 3 , 20.8% 
for tetragonal Zr 02 . Notice that these values are close 
to, but do not exactly match, those entering the defini¬ 
tion of standard hybrid functionals, such as PBEO (25%) 
and B3LYP (20%) Also, the corresponding dielec¬ 
tric constants are found close to the experimental values, 
deviating from them by at most four percent (see also 
Refs. [H andiiol). 

C. Charge-transition levels 

The position of the electronic levels introduced by O 
vacancies in the band gap of the material was computed 
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on the basis of the CTLs formalism^"— The latter con¬ 
stitutes a rigorous framework for computing excitation 
and emission energies in defective semiconductors that 
can be directly compared with experiments (see the In¬ 
troduction) . 

At fixed experimental conditions, corresponding for in¬ 
stance to exposing the material to high-temperature O- 
rich or 0-poor atmospheres in the case of oxide-based 
materials, the formation energy of the defect (O vacancy) 
depends only on the chemical potential (which at T = 0 
and without an external field can be identified with the 
Fermi energy) of the electrons that can be exchanged 
between the defect and the host crystal subsystems, fol¬ 
lowing some external (e.g., optical) perturbation. 

The optical transition level is defined as the 

Fermi energy at which the formation energies of the de¬ 
fect in the (dimensionless) charge states q, and g', 

Aq , become equal. Since we are dealing with opti¬ 
cal (vertical) excitations, occurring over time-scales much 
shorter than the phonon time-scale, the atomic configu¬ 
ration is kept frozen at that of the initial defect charge 
state, and all the quantities of interest (e.g., the forma¬ 
tion energy) are evaluated at that configuration for both 
charge statesi^ The optical transition level, referred to 
the top of the valence band (VB), £v, of the host crystal, 
can thus be expressed as 

( 1 ) 

q-q 

where the the ground-state total energies of the defective 
supercells with defect in the charge states q and q', Eu,q 
and Ed^, are presently assumed to incorporate any cor¬ 
rections due to electrostatic interactions between charged 
image defects and to the definition of a common energy 
reference for the two bulk calculations 

Using Janak’s theorem,— and in the spirit of Slater’s 
transition-state theory^ the total energy difference ap¬ 
pearing in Eq. © can be approximated as^^i^ 

EB,q' - Eu.q ^ ^ [sn{N; q') + eN{N -l;q = q +1)], 

( 2 ) 

where £NiN; q') is the KS eigenvalue of the highest occu¬ 
pied spin-orbital of the A^-electron system (corresponding 
to defect charge state q') and e n{N — 1] q = q' -\-1) is the 
KS eigenvalue of the lowest unoccupied spin-orbital of 
the {N — l)-electron system, corresponding, for instance, 
to excitation of an excess electron into the host crystal 
(defect charge state q = q' 

Two issues have to be considered in the evaluation of 
the KS eigenvalues in Eq. ©: 

• The KS eigenvalues obtained from the two bulk cal¬ 
culations (defective supercells corresponding to dif¬ 
ferent charge states q and q') have to be referenced 
to some common energy scale. We chose the zero of 
the energy to be the Is KS eigenvalue of O (located 
~ 500 eV below the VB edge), averaged over the 
whole supercell. 


• For charged defects, the electrostatic interactions 
between image defects arising from usage of finite- 
size supercells have to be taken into account.— 
The corresponding correction to the KS defect level 
Ed reads^ 

AeD(Ar;g) = -^A<P = -2(l + /)(|^)^, (3) 

where the correction AU^^ to the total energy of 
the charged supercell is evaluated in the Makov- 
Payne schemej^i^ om is the Madelung constant, L 
is the average distance between image defects, and 
Coo is the static electronic dielectric constant of the 
host crystal, computed within the self-consistent 
hybrid scheme— The prefactor / Ri —0.35 approx¬ 
imately takes into account the interaction energy 
of the defect charge with the neutralizing jellium 
background (scaling as L~^) introduced to renor¬ 
malize total energy in charged defect calculations— 

Thermodynamic levels /r*’’'®™(g/g'), corresponding to 
excitations qt —^ qf in which the system is allowed to 
relax in the final charge state qf (adiabatic transitions) 
are obtained from optical levels by adding or subtracting 
the relaxation energyi^ 

Erel = 1^. - Eu,qj , (4) 

according to whether q = qf = qi + 1 = q' + 1 (elec¬ 
tron excitation) oi q' = qf = qi — 1 = q — 1 (electron 
trapping)— Notice that in Eq. ([4]) electrostatic correc¬ 
tions to formation energies would have to be evaluated 
for the same charge state g/, and they thus cancel against 
each other- 


ill. RESULTS AND DISCUSSION 
A. Ti02 

The optoelectronic properties of 0-deficient Ti02 
(Ti 02 -a;) have been extensively investigated both 
experimentally^ and theoreticallyr^i^^— making this 
system an ideal model for testing the performance of 
the dielectric-dependent hybrid method in describing the 
electronic properties of defective semiconductors. 

The observed n-type conductivity of Ti02 samples has 
been related to the presence of O vacancies, acting as 
intrinsic donors, formed under strongly reducing synthe¬ 
sis conditions— The larger intrinsic conductivity of the 
anatase phase, which also makes it more appealing for 
applications, has been explained with the different be¬ 
havior of the O vacancy in the two polymorphs— 

On the other hand, optical experiments found the as¬ 
sociated defect levels to be located deep in the band gap. 
Infrared absorption (IR)j^ ultraviolet photoelectron^ 
and electron energy-loss experiments^ on reduced rutile 
samples all agree in attributing a feature at about 1 eV 
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below the conduction band (CB) to the presence of O va¬ 
cancies. Theoretical calculations have revealed that the 
apparent contradiction between the observed n-type con¬ 
ductivity and the deep nature of optical levels in Ti 02 -a; 
is rationalized by the different character (shallow versus 
deep) of the 0-vacancy thermodynamic and optical tran¬ 
sition levels.— 

Furthermore, the excess charge due to O deficiency has 
been shown to redistribute differently in the host crystal, 
according to whether the rutile or anatase polymorph 
is considered, resulting in the localization of the extra 
electrons on topologically different Ti atoms— 

It is the goal of the present part of the work to re¬ 
investigate these issues within the dielectric-dependent 
hybrid approach, which by construction is sensitive to the 
different electronic screening characteristic of the rutile 
and anatase phases. 


1. Rutile 


The ground state of rutile Ti02 with a neutral O va¬ 
cancy is found to be a triplet, with the two excess elec¬ 
trons redistributing over the whole (110) plane contain¬ 
ing the vacancy [Figure [IJa)]. Most of the excess charge 
is contributed by the fully-coordinated Ti 2 (spin den¬ 
sity Us = = 0.70), Tia (two equivalent atoms. 

Us = 0.15 each) and T^ {us = 0.65) atoms. The 
two under-coordinated Tii atoms instead accommodate 
a smaller fraction of it {ug = 0.13 each). This picture 
is in qualitative agreement with recent scanning tunnel¬ 
ing microscopy (STM) measurements on the rutile (110) 
surface^ suggesting the localization of the extra elec¬ 
trons on Ti atoms away from the vacancy^ with the 
formation of a polaron. Electron paramagnetic resonance 
(EPR) measurements also confirmed the presence of Ti^’*' 
ionic species formed as a consequence of electron trapping 
at Ti'^^ centers in Ti 02 -x crystals— 

Notice that the localization of the unpaired electrons 
at Ti^+ centers is correctly described only within the 
hybrid-functional scheme;^ at variance with the case 
of standard local/semilocal functionals which commonly 
fail, due to incomplete cancellation of the self-interaction 
error. Notice also that the above physical picture is cor¬ 
rectly captured only at the level of spin-polarized DFT, 
while neglecting spin polarization leads to a completely 
different description, with the excess charge trapped in 
the vacancy void )^^'^^ in a configuration typical of a color 
center (F center). We found the corresponding closed- 
shell ground state to be highly unstable, being above 
the triplet solution by nearly 2 eV in total energy. In 
the spin-compensated case, the two under-coordinated 
Ti atoms relax towards the vacancy, shortening their dis¬ 
tance to 2.90 A (compared with an unrelaxed distance of 
2.99 A); instead, in the true, spin-polarized ground state 
(triplet), they move away from the vacancy, resulting in 
an equilibrium distance of 3.37 A. 

Our found magnetic ground state is compatible with 



FIG. 1. O vacancy in rutile Ti02. (a) Isosurface plot of the 
charge density associated with the excess electrons for the 
neutral vacancy. The yellow (large) and red (small) spheres 
represent Ti and O atoms, respectively. Only the atomic 
plane containing the O vacancy is shown, (b) Correspond¬ 
ing electronic band structure (only spin majority is shown) 
computed at the equilibrium geometry, (c) Optical (solid, 
black) and thermodynamic (dashed, red) charge-transition 
levels. Positions (in eV) are given with respect to the top 
of the valence band (VB). 


the evidence of triplet and doublet signals obtained in 
EPR measurements^2iS£ on rutile Ti02-a: crystals, and 
attributed to neutral and singly-charged O vacancies, re¬ 
spectively. 

We succeeded in stabilizing a triplet configuration in 
which nonetheless the excess charge is confined at the 
vacancy void;^ a situation typically encountered in non¬ 
reducible 0-deficient oxides— However, we found the 
corresponding ground-state energy to be ^ 0.8 eV higher 
than that of the triplet polaron-like solution represented 
in FiguredKa). This finding definitely rules out the possi¬ 
bility of an F-center-like behavior of O vacancies in rutile, 
in favor of the polaron picture, which is in agreement with 
all the available experimental evidence. This conclusion 
is also corroborated by previous theoretical investigations 
using larger supercells— 

The bulk band gap of rutile Ti02, as computed us¬ 
ing the self-consistent hybrid functional, is found to be 
2.99 eV, a value compatible with the available data from 
photoemission/inverse photoemission experiments, yield¬ 
ing an electronic gap of 3.3±0.5 eV (Ref.|6l) and ~ 3.1 eV 
(Ref. [6^. Upon removal of an O atom, the two resulting 
excess electrons occupy two triplet defect states appear- 
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ing in the band gap at about 0.5 — 1 eV below the CB 
[Figure [TJb)]. As the vacancy is ionized, the remaining 
excess electron localizes on the three Ti atoms in the row 
of the O vacancy, again with a dominant contribution 
from the fully-coordinated Ti 2 atom (ris = 0.68). The 
associated occupied defect level is still localized at about 
1 eV below the CB. 

Figure [T][c) shows the computed charge-transition lev¬ 
els. The optical levels lie deep in the band gap, resulting 
in excitation energies from the defect states to the CB of 
0.47 eV, 1.20 eV and 0.84 eV for the (-hl/0), (+2/ -h 1) 
and (-1-2/0) transitions, respectively. Notice that, as a re¬ 
sult of the finite supercell size, the residual defect-defect 
interaction (i.e. overlap between polarons)^ leads to 
a small dispersion of the highest-occupied defect band 
(about 0.2 eV). This should be considered as indicative 
of the numerical accuracy of the computed level. Be¬ 
ing aware of this limitation, our results can be consid¬ 
ered consistent with the features, found at about 1 eV 
below the bottom of the CB, observed in experimental 
spectrai^^— In particular, IR measurements^ revealed 
two absorption peaks at 0.75 eV and 1.18 eV for the sin¬ 
gle and double ionization, which correlate well with our 
computed first and second optical ionization energies for 
the O vacancyi^ 

Thermodynamic transition levels are significantly 
higher in energy, reflecting the sizable structural re¬ 
laxation associated with the polaronic distortion. The 
(-1-1/0) adiabatic transition is located about 0.1 eV be¬ 
low the CB, indicating that the neutral O vacancy may be 
stable at sufficiently low temperature. The second tran¬ 
sition (-1-2/-1-1) is at 0.81 eV, while the double ionization 
[(-1-2/0) transition] requires 0.46 eV. These values should 
be compared with measured thermal ionization energies 
observed in a range between 0.3 and 0.6 eVi^ 



FIG. 2. O vacancy in anatase Ti02. (a) Isosurface plot of the 
charge density associated with the excess electrons for the 
neutral vacancy. The yellow (large) and red (small) spheres 
represent Ti and O atoms, respectively, (b) Corresponding 
electronic band structure (only spin majority is shown) com¬ 
puted at the equilibrium geometry, (c) Optical (solid, black) 
and thermodynamic (dashed, red) charge-transition levels. 
Positions (in eV) are given with respect to the top of the 
VB. 


2. Anatase 

The nature of the neutral O vacancy in anatase 
Ti 02 has been extensively studied in the theoretical 
literature42i^^i^Sii However, no consensus has been 
achieved so far concerning the redistribution of the asso¬ 
ciated excess electrons in the host crystal. This difficulty 
may be rationalized with the presence of several minima 
on the relevant adiabatic potential energy surface. The 
realization of a specific configuration is strongly related 
to the description of the structural relaxation around the 
vacancy, which is affected, for instance, by the employed 
DFT xc approximation. 

Our calculations yield a triplet ground state for the 
system with a neutral O vacancy, with a correspond¬ 
ing substantial rearrangement of the atomic positions 
around the vacancy. As shown in Figure [2][a), this im¬ 
plies a reduction of the local symmetry around the va¬ 
cancy, with a nearest-neighbor O atom breaking its bond 
with a bulk Ti atom (Ti 2 ) and considerably relaxing to¬ 
wards the vacancy void^^^ The same behavior was ob¬ 


served in all previous calculations in which symmetry¬ 
breaking relaxation was allowed fori ^s,5i,69 qijjg 
cess electrons are found to be localized on two Ti atoms 
(the under-coordinated Tii, with Ug = 0.92, and the 
bulk Ti 2 , with Ug = 0.84), similarly to a previous 
DFT-I-U investigation!^ This finding is again in quali¬ 
tative agreement with the STM investigation of Diebold 
and coworkers,— which suggests that the excess charge 
preferably stays close to the vacancy in anatase, as op¬ 
posed to the case of rutile in which the major contri¬ 
bution is given by the fully-coordinated Ti atoms (see 
Section lill A II) . 

For the system with a neutral vacancy, the computed 
electronic structure shows two triplet defect states at 
^ 0.7—1 eV below the CB [Figure [2][b)]. The com¬ 
puted band gap is 3.69 eV, larger than the value of 
3.42 eV obtained from optical absorption experiments 
at low temperature— To the best of the authors’ knowl¬ 
edge, no data for the photoemission gap of stoichiometric 
anatase Ti 02 are available in the experimental literature. 
However, our computed gap is consistent with previous 
GW calculations— 
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When the vacancy is ionized, a doublet occupied de¬ 
fect state locates at ^ 0.7 eV below the CB. The ex¬ 
tra electron remains localized (rig = 0.92) at the under¬ 
coordinated Til atom. 

The computed thermodynamic transition levels are sig¬ 
nificantly shallower than in rutile [compare Figure [ijc) 
and Figure [2l[c)]. In particular, the (-1-1/0) level is prac¬ 
tically resonant with the CB, suggesting that an electron 
can be thermally excited into the CB, thus accounting 
for the observed larger conductivity of anatase samples. 

Similarly to the case of rutile, optical levels are 
substantially deeper than the corresponding thermo¬ 
dynamic ones. We computed excitation energies of 
0.64 eV, 1.15 eV and 0.89 eV for the (-fl/O), {+2/ + 
1) and (-1-2/0) transitions. Experimentally, a feature 
at ^ 1.0 — 1.1 eV below the CB has been reported 
in different experiments (combined x-ray photoemis¬ 
sion/absorption spectroscopy;^ as well as scanning tun¬ 
neling spectroscopy^) on reduced anatase surfaces, and 
has been attributed to O vacancies. Notice that, since the 
stable charge state of the O vacancy is g = -1-1 at room 
temperature, it seems likely that the (- 1 - 2 / -|- 1 ) transition 
is mainly probed in optical experiments; accordingly, our 
computed excitation energy of 1.15 eV matches excel¬ 
lently with the measured spectroscopic features. 

In conclusion, the present approach proved capable of 
explaining several experimental signatures related to the 
different behavior of O vacancies in the two most common 
Ti 02 polymorphs, rutile and anatase. In particular, the 
larger conductivity of anatase can be understood by anal¬ 
ysis of the thermodynamic transition levels, whereas the 
deep nature of the electronic features observed in spec¬ 
troscopies is rationalized on the basis of the computed 
optical transition levels. Our findings are in qualitative, 
and sometimes quantitative, agreement with previous in¬ 
vestigations at the hybrid- functionali^i^ and DFT-|-U^ 
level. 


B. WO3 

The optical and electrical properties of substoichio- 
metric WO 3 (WOa-x) have been extensively investi¬ 
gated experimentally, particularly in relationship with 
the electrochromism exhibited by this materiali^ How¬ 
ever, the microscopic origin of this effect, which is in¬ 
timately connected with the typical blue coloration of 
WOs-a; filmshas been longly debated, and several 
models have been put forward for rationalizing experi¬ 
mental observations 

On the theory side, a few first-principle investiga¬ 
tions have been performed on 0-deficient RT monoclinic 
WOa;^^— but among them only one^ employed state- 
of-the-art methods yielding a correct description of the 
bulk electronic structure. The latter study of Wang et 
al^ revealed a delicate interplay between the concentra¬ 
tion of O vacancies and the metallic or insulating nature 
of WOs-a;. Furthermore, different behaviors were ob- 
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FIG. 3. O vacancy in RT monoclinic WO 3 . Isosurface plot 
of the charge density associated with the excess electrons for 
the (a) neutral and (c) singly-charged vacancy. The yellow 
(large) and red (small) spheres represent W and O atoms, 
respectively. Corresponding electronic band structures, com¬ 
puted at the equilibrium geometry, for the (b) neutral and (d) 
singly-charged vacancy (only spin majority is shown). (e) Op¬ 
tical (solid, black) and thermodynamic (dashed, red) charge- 
transition levels. Positions (in eV) are given with respect to 
the top of the VB. 


served according to the orientation of the W-O-W chain 
along which the O atom is removed. For the sake of sim¬ 
plicity, we will limit the present discussion to the case of 
an O vacancy created along the a crystallographic axis 
of the monoclinic cell;^ which we have evidence to be¬ 
lieve to be the main responsible for the observed optical 
features of WOs-a,, namely its blue coloration. A more 
comprehensive investigation will be subject of a future 
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publication. 

Our calculations indicate the closed-shell singlet solu¬ 
tion to be more stable (by 0.17 eV) than the open-shell 
triplet one for the system with a neutral vacancy. Corre¬ 
spondingly, the two extra electrons occupy a defect state 
located at about 1.0 —1.5 eV below the CB [Figure[3Ka)]; 
they localize mainly at the vacancy void, with a contri¬ 
bution from the two under-coordinated W atoms [Fig¬ 
ure [SKb)]. 

As the vacancy is ionized, the extra electron localizes 
on one of the under-coordinated W atoms [Figure [21[c)], 
leading to the formation of a reduced ionic state. In 
fact, one of the proposed models for the chromic mecha¬ 
nism in WO 3 -X suggests that the presence of cen¬ 
ters results from the localization of excess electrons on 
sites upon light absorption^S As shown in Fig¬ 
ure |3][d), this 5d electron occupy a defect state ^ 1.3 eV 
below the CB 2 eV if the geometry of the neutral 
vacancy is retained). 

Several experiments have been performed to elucidate 
the spectroscopic properties of substoichiometric WO 3 
samples. The coloration efficiency of WOs-a, {x < 0.4) 
films has been reported to increase with increasing O 
deficiency;^ Their typical blue coloration has been cor¬ 
related with a broad absorption band with a maximum 
at 900 nm 1.38 eV)^!^^ although this feature has 
been observed in amorphous films. Photoluminiscence 
(PL) spectra^^ exhibited an emission peak at 550 nm 
(^ 2.26 eV), and photoelectron measurements22i2^ con¬ 
firmed the presence of a defect state attributed to O va¬ 
cancies at ^ 2 eV above the VB. 

The computed optical transition levels reported in Fig¬ 
ure He) are consistent with these observations, being po¬ 
sitioned in the band gap in the range from 1.72 eV to 
2.33 eV above the VB. The dispersion of the neutral- 
vacancy defect state limits the accuracy of the calculated 
transition levels to ~ 0.2 eV. Notice that the computed 
band gap is larger than the reported values from absorp¬ 
tion measurements (see Ref. l2J and references therein), 
but in excellent agreement with ultraviolet direct/inverse 
photoemission experiments (3.45 eV from Ref. [oO), as 
well as with GW calculationsThus, the excitation 
energy pertaining to, e.g., the (- 1 - 1 / 0 ) transition, which 
is obtained to be 1.12 eV, is compatible with the experi¬ 
mentally observed features in the absorption spectra. 

The corresponding thermodynamic levels are predicted 
to be well-detached from the bottom of the CB, with a 
minimum excitation energy of 0.46 eV computed for the 
(-1-1/0) transition. This implies that the specific kind of 
O vacancy considered here cannot contribute to the ob¬ 
served n-type conductivity of WOa-x films, which has 
been attributed to O vacancies;^ Thus, although the 
present work clarifies the nature of the spectroscopic fea¬ 
tures commonly evidenced in reduced WO 3 , in agreement 
with the previous investigation of Wang et al.^ further 
study is called for to elucidate the possibly different be¬ 
havior of other inequivalent O vacancy sites. 


TABLE I. Formation energy for the O vacancy in the different 
materials, computed with respect to ( 1 / 2 )A[ 02 ], where i?[ 02 ] 
is the ground-state total energy of the triplet O 2 molecule. 



X 

Formation energy (eV) 

Ti 02 -i (rutile) 

1/24 

5.2 

Ti 02 -i (anatase) 

1/32 

4.7 

WO3-. 

1/16 

5.3 

Z1O2-X 

1/36 

6.6 


C. Zr02 

The O vacancy in reduced zirconia (Zr 02 -a;) has been 
reported to exhibit features similar to those of the F 
centers typically observed in non-reducible oxides^ the 
extra charge being trapped in the vacancy voidi^^ EPR 
studies^^ indicate that Zr 02 is generally less prone to 
lose oxygen than other reducible oxides, such as Ti 02 . 
This is confirmed by our calculations, giving a formation 
energy for the O vacancy in Zr 02 -a; larger by more than 
1 eV with respect to the other oxides (see Table © , in 
agreement with previous investigations 

As shown in Figure 0Ka), upon removal of an O atom 
in Zr 02 , the two excess electrons stabilize in the vacancy 
void, with a contribution from the Ad orbitals of two 
nearest-neighbor Zr atoms (Zri and Zr 2 ); the correspond¬ 
ing ground state is a closed-shell singlet. As expected, the 
neighboring Zr atoms relax towards the vacancy, result¬ 
ing in a shrinking of the distance Zri-Zr 2 from 3.61 A 
to 3.48 A, and of Zri-Zr 3 from 3.65 A to 3.61 A. The 
corresponding band structure [Figure|ll[b)] exhibits a flat 
defect state positioned at ^ 3.3 eV above the VB, in close 
agreement with a previous B3LYP studyi^ 

As the vacancy is ionized, the remaining extra electron 
distributes on the two Zri and Zr 2 atoms, as well as in the 
space between them, in a typical bonding configuration. 
The extra electron occupies a doublet state in the band 
gap at ^ 3.6 eV above the VB. 

We also studied the case of the negatively-charged O 
vacancy, which can be realized as a consequence of the 
trapping of an electron from the bulk CB. Such trapping 
phenomena can affect, for instance, the performance of 
electronic devices, for which high-K oxide materials such 
as zirconia or hafnia (Hf 02 ) have been proposed as gate 
electrodesiSi The trapped electron occupies an additional 
defect state in the band gap at ~ 1.2 eV above the (now 
spin-resolved, with spin majority shown) states associ¬ 
ated with the neutral vacancy [Figure|U[d)], in a bonding 
configuration between the Zri and the Zr 2 atoms [Fig¬ 
ure (He)]. 

The computed CTLs corresponding to both vacancy 
ionization and electron trapping are shown in FigureSKe) . 
A band gap of 5.73 eV is calculated, to be compared with 
the experimental value of 5.78 eV deduced from optical 
absorption measurements;^ typically GW calculations 
give a band gap larger by ^ 0.1 — 0.2 eV (see Ref. m and 
references therein). 








FIG. 4. O vacancy in tetragonal Zr02. Isosurface plot of the charge density associated with (a) the two excess electrons for the 
neutral vacancy, and (c) the trapped electron for the negatively-charged vacancy. The yellow (large) and red (small) spheres 
represent Zr and O atoms, respectively. Corresponding electronic band structures, computed at the equilibrium geometry, 
for (b) the neutral, and (d) the negatively-charged vacancy (only spin majority is shown), (e) Optical (solid, black) and 
thermodynamic (dashed, red) charge-transition levels. Positions (in eV) are given with respect to the top of the VB. 


Much experimental work is available in the literature 
investigating the PL properties of Zr02-a;. However, a 
clear understanding of the underlying mechanisms seems 
to be lacking. Several emission features in the region be¬ 
tween 2.0 and 3.5 eV are typically observed, and have 
been related to the presence of O vacancies)^ The large 
uncertainty in experimental results may be the conse¬ 
quence of the practical difficulty in obtaining high-purity 
single crystals of zirconia with a well-defined phase. 

PL measurements on tetragonal nanocrystalline 
Zr02-x showed an emission peak at 350 nm (^ 3.54 eV), 
which was attributed to electron-hole recombination in¬ 
volving F-center states in the band gap. Analysis of the 
PL spectrum evolution with annealing allowed to assign 
this feature to O vacanciesi^ Our computed position of 
the (-1-1/0) optical transition level (3.40 eV above the top 
of the VB) may account for this observation. 

Another study on undoped tetragonal zirconia 
nanocrystals^^ evidenced the presence of an emission 
peak at ~ 2.8 eV. This may be related to a PL process 
involving the (-1-2/ -|- 1) optical level, which we found at 


2.77 eV above the top of the VB. 


The thermodynamic levels are found below the CB by 
at least 1.5 eV, indicating that all the 0-vacancy charged 
states are stable at room temperature. An exception 
is represented by the (0/ — 1) level, which is separated 
from the CB by ~ 0.3 eV. This level is associated to 
electron trapping processes causing the 0-vacancy to be¬ 
come negatively-charged. Our calculations suggest that 
such trapping center is stable at room temperature, and 
the same conclusion was obtained from a similar investi¬ 
gation on hafniai^^ in which, due to analogous chemical 
behavior of the Hf and Zr cations, O vacancies are ex¬ 
pected to induce similar effects as in zirconiajiSS In this 
respect, trapping/detrapping experiments on Hf 02 found 
the activation energy for the trap level to be 0.35 eVfiSi 
which is in agreement with our calculations. 
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IV. CONCLUSIONS 

We reviewed the behavior of the O vacancy in proto¬ 
typical reducible (Ti 02 , WO 3 ) and non-reducible (Zr 02 ) 
oxide semiconductors adopting a novel DFT computa¬ 
tional approach based on hybrid functionals with the 
exact-exchange fraction obtained ab initio for the bulk 
non-defective material. Its capability of yielding accu¬ 
rate electronic band structures, together with structural 
properties and total energies, makes it an ideal tool for 
investigating defect levels in semiconductors. 

The computed optical transition levels are generally 
found in agreement with various spectroscopy experi¬ 
ments, providing support for their interpretation. In¬ 
stead, thermodynamic levels allow one to infer about the 
intrinsic degree of n-type conductivity of substoichiomet- 
ric oxide materials through analysis of the stability of the 
different 0-vacancy charge states. Residual discrepancies 
between theory and experiment may arise from two main 
sources: (i) the always present approximate treatment of 
xc terms, and (ii) usage of finite supercells. While the 
former is intrinsically related to the adopted functional, 
the latter may show up also as consequence of the evalua¬ 
tion of CTLs via the Janak’s theorem, Eq. (ED, instead of 
by direct computation of total energy differencesdue 
to the defect band dispersion introduced by the residual 
polaron overlapping in neighboring cells. 

Fundamental differences are evidenced in the behav¬ 
ior of reducible and non-reducible oxides. Upon creation 
of an O vacancy, the former tend to exhibit polaron-like 
features: the excess charge mainly localizes on the metal¬ 
lic ions, reducing their formal charge, and the relaxation 
of the atomic environment is not restricted to the atoms 
in the immediate vicinity of the defect. Instead, in non¬ 
reducible oxides such as Zr 02 , the excess charge is sta¬ 
bilized in the vacancy void, and the relaxation is more 
local; furthermore, transition levels are deeper, being sit¬ 
uated near the center of the band gap. These are the 
fingerprints of the F-center-like behavior of O vacancies 
in these materials. 

In conclusion, we believe that the present approach 
is able to characterize defect levels in semiconductors 
in a way that makes it valuable for predictive stud¬ 
ies also on extrinsically-doped semiconductors. Notice 
however that, while the method is straightforwardly ap¬ 
plied to point (zero-dimensional) defects, its extension 
to the treatment of higher-dimensional defect structures, 
such as surfaces, requires facing unprecedented difficul¬ 
ties. For instance, modeling the interface between two 
different semiconductors implies the exchange fraction 
to become a space-dependent quantity, and this even¬ 
tually calls for a phenomenological model of the optical 
properties of the interface. To the best of the author’s 
knowledge, the issue has not yet been considered in the 
literature, and its solution would make it possible to ap¬ 
ply the method to investigate the rich physical-chemical 
phenomena occurring at interfaces. 
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